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The Landau gauge lattice gluon propagator is discussed for different sets of lattices. Particular at- 
tention is given to its infrared properties. Our results show that the lattice propagator can be made 
compatible with either the decoupling-like or the scaling-like solution of the Dyson-Schwinger 
equations. Furthermore, the analysis of the Cucchieri-Mendes bounds is performed considering 
large volume simulations and the Oliveira-Silva ratios are computed. If the first do not give a 
clear answer about the value of D(Q), the second method favors a D(Q) = 0. Finally, the SU(3) 
and SU(2) propagators are compared in the infrared. It comes out that the propagators are differ- 
ent although the infrared exponents seem to be similar. The analysis suggests a scaling behaviour 
D(0) ~ N with the gauge group SU(N). 
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1. Introduction 

For a given quantum field theory the Green's functions encode the dynamical information. For 
QCD, in particular, the computation of any Green's function, such as the gluon propagator D(q 2 ), 
over the entire momentum spectrum cannot rely on perturbation theory. In what concerns the gluon 
propagator, in the past years, there has been a discussion about its infrared behaviour and its value 
at zero momentum in the Landau gauge. The recent effort on computing D(0) comes from its 
relation with the Gribov-Zwanziger gluon confinement mechanism, which implies D(0) = 0, and 
with the possibility of dynamical mass generation for gluons. The dispute is still going on and 
involves both Schwinger-Dyson solutions and lattice QCD results. 

The recent solutions of the Dyson-Schwinger equations have different infrared behaviours. On 
one side we have the scaling solution |Jlj, ||] with D(0) = and an infrared behaviour given by a pure 
power lawD(4 2 ) = (q 2 ) 2K -\ with K ~ 0.595. On the other side we have the decoupling solution 
|j] (see also ^) with a finite and nonvanishing D(0), with the value of D(0) being related with 
a dynamical generated gluon mass, and a plateau for D{q 2 ) at low momenta. 

The recent lattice simulations also show contradictory results. Indeed, large volume simula- 
tions using the Wilson action and large lattice spacings, i.e. a ~ 0.18 fm or larger, show a gluon 
propagator that agrees qualitatively with the decoupling solution. However, as we will see, a sec- 
ond look at the lattice data seems to indicate that the solution D(0) =0 is still compatible with the 
lattice data — see below. 

Here, we want to discuss on what the lattice simulations tell us about the infrared gluon prop- 
agator and, hopefully, point towards a "favorite" value of 0(0). In the following, we will use 
standard definitions for each of the quantities, which will not be shown here. The interested reader 
can find the details for example in [§]. 

2. Lattice setup 

Our discussion of the infrared gluon propagator in the Landau gauge uses lattice data for the 
Si/ (3) gluon propagator from simulations with two different values of j8, namely j8 = 6.0 and 
j8 =5.7, combining data generated at Coimbra with the data from the Berlin-Moscow- Adelaide 
group Furthermore, the SU(3) data will be compared with the large volumes, i.e. 80 4 and 128 4 
at /3 = 2.2, SU(2) propagator of the S. Carlos group [[J. 

The lattices simulated using the SU(3) gauge group and the number of configurations for 
p = 5.7 are 

L 8 10 14 18 26 36 44 64* 72* 80* 88* 96* 
L(fm) 1.5 1.8 2.6 3.3 4.8 6.6 8.1 11.8 13.2 14.7 16.2 17.6 
#Confs 56 149 149 149 132 100 29 14 20 25 68 67 

where * stands for simulations carried out by the Berlin-Moscow- Adelaide group Note that we 
took the lattice spacing from the string tension |7|] and not from ro as the Berlin-Moscow-Adelaide 
group did in [JSj] . We have rescaled their data according to our definitions. We would like to call the 
reader's attention that different volumes have different statistics. The lattice gluon propagator for 
j3 = 5.7 can be seen in figure [|. 
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Figure 1: Bare gluon propagator from j3 = 5.7 simulations. The data for the largest volumes at L > 64 is 
the Berlin-Moscow-Adelaide data [Kfl rescaled to our definition of the lattice spacing. 
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Figure 2: Bare gluon propagator from j3 = 6.0 simulations. 



The simulations for j8 = 6.0, where the lattice spacing is a = 0.1016(25) fm, were carried for 
the following lattices 

and the propagators are plotted in figure |[ 

As seen in figures [I] and 0, no supression of the gluon propagator is observed. Therefore, as a 
first and naive conclusion one could claim that, modulo finite volume effects, the zero momentum 
gluon propagator is finite and non-zero. 



3 



Lattice Infrared Landau Gauge Gluon Propagator 



O. Oliveira 



L 16 20 24 28 32 48 64 80 
L(fm) 1.6 2.0 2.4 2.8 3.2 4.9 6.5 8.1 
#Confs 52 72 60 56 126 104 120 18 



3. Modelling The Gluon Propagator 

As a first step towards trying to distinguish between a vanishing or nonvanishing D(0), we look 
at the compatibility of the lattice data with the functional forms which have been used to describe 
the two Dyson-Schwinger solutions, i.e. the infrared lattice data is fitted to 

where M 2 plays the role of a hard mass, and 

D{q 2 ) = Z-i±J (3.2) 

(q 2 + A 2 ) 2K 



In the fits to ( p.2j ) the point q = is not included. This is not unwise since finite volume effects are 
certainly larger for 0(0). 

In what concerns the fits to the /3 = 5.7 data, we have observed that for the massive like 



propagator (3.1), the smaller lattices (L < 6fm) are not described by the above functional form. 
Moreover, for the largest lattices, (3.1) reproduces well the lattice data, i.e. fits have % /d.o.f. < 



1.8, for momenta up to 500 MeV with a gluon mass in the range of 719(18) MeV to 2.20(63) GeV. 
M depends on the fitting range and on the lattice volume. It was observed that M (and Z) tend to 



increase with the lattice volume. For the fits to ( |3.2[ ), the results are similar. Indeed, the lattice 
data is well described by ([0|) to momenta up to 500 MeV, with A decreasing from ~ 900 MeV 
to ~ 420 MeV as the lattice volume increases. The K is, within one standard deviation, above 0.5 
suggesting that D(0) = 0. Typical fits can be seen in figures |3| and ||. 

We repeated the procedure using the SU(2) gluon propagator of [^] and found a similar be- 
haviour. The main difference being that the SU(2) mass scales are typically larger than the corre- 
sponding SU(3) mass scales. For M we get numbers around 1 GeV, while A stays between 700 
MeV to 800 MeV and K touch values slightly below the "magics" 0.5. 

The results of fitting the j8 = 6.0 lattice gluon data give similar results. For completeness we 
resume the results as follows: 

1. "Decoupling" Fit: 

48 4 (L=4.9 fm), q up to 508 MeV; M goes from 600(20) MeV down to 578(12) MeV as the 
fitting range is increased; 

64 4 (L=6.5 fm), q up to 503 MeV; M goes from 753(57) MeV down to 655.8(9.3) MeV as 
the fitting range is increased; 

80 4 (L=8.1 fm), q up to 664 MeV; M goes from 588(129) MeV down to 576.4(6.6) MeV as 
the fitting range is increased; 
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IR Bare Gluon Propagator - p = 5.7 
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Figure 3: IR bare gluon propagator for /3 = 5.7 and fits to (2.1). 
IR Bare Gluon Propagator - p = 5.7 
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Figure 4: IR bare gluon propagator for j3 = 5.7 and fits to (2.2). The MeV scales in the graph show the 
maximum fitting range compatible, i.e. with a X 2 /d.o.f. < 1.8, with (2.2). Note that for the largest volume 
the q = GeV point is missing. 

2. "Scaling" Fit: 

48 4 (L=4.9 fm), q up to 671 MeV; A goes from 460(61) MeV down to 432(27) MeV and k 
goes from 0.579(56) to 0.606(34) as the fitting range is increased; 

64 4 (L=6.5 fm), q up to 503 MeV; A goes from 609(111) MeV up to 614(43) MeV and k 
goes from 0.510(31) to 0.513(16) as the fitting range is increased; 

80 4 (L=8.1 fm), q up to 645 MeV; A goes from 1.1(1.8) GeV down to 582(54) MeV and k 
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IR Bare Gluon Propagator - p = 5.7 
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Figure 5: Exponent K as measured from fitting (2.2) (black points). The plot also shows K measured using 
/3 = 6.0 asymmetric lattice data to fit (2.2) (blue line) and ratios of propagators (green line) using the same 
asymmetric lattice set. The red line is K = 0.5 and a value of K above 0.5 implies D(Q) = 0, a K < 0.5 means 
D(Q) = °o and K = 0.5 implies a decoupling type solution. 



goes from 0.453(61) to 0.556(17) as the fitting range is increased l . 

From the previous analysis, the conclusion is that although a "decoupling" like gluon propaga- 
tor seems to be favoured by the raw lattice data, a "scaling" like propagator is not excluded yet. The 
fits just described show that both type of solutions are in good agreement with the lattice data and 
that if a "scaling" like propagator is the solution, D{q 2 ) starts to be supressed only for rather low 
momenta. Believing on the results of figure D{q 2 ) is supressed only for momenta well below 
100 MeV. 



4. Cucchieri-Mendes Bounds 

In JTo|], the authors derived inequalities between D(0)/V and what they called an average 
absolute value of the components of the colour magnetization M(0), 

(M(0)) 2 < < d {N 2 C - 1) (M(0) 2 } , (4.1) 

where d is the number of space-time dimensions and N c the number of colours. For the definition 
of M(0) see the cited work. In the above expression, (• • •) means Monte Carlo average and ( fi~Tl ) 
follows directly from using a Monte Carlo approach. 



In [10] the different terms in (4.1 ) were computed for SU(2) simulations and, after performing 



a scaling analysis, the authors conclude in favour of a finite and nonvanishing D(0). Indeed, they 



'The reader should be aware that the presented results for the 80 4 lattice have a quite small statistics. So far, the 
total number of gauge configurations is 18, which could explain the relative large errors in A and the rather low value 
for k obtained with the smallest fitting range. 
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Figure 6: Cucchieri-Mendes bounds, the lattice data and corresponding fits to A/V a . Note that in all fits 
the 26 4 data were excluded in order to have a % 2 I d.o.f . < 1 .8. Furthermore, for D(0) /V the plot includes 
the fits to the smaller set of lattices, i.e. 8 4 — 44 4 and to all data available, i.e. to the lattices 8 4 — 88 4 . The 
agreement between the two sets of data is perfect. 



claimed D(0) > 2.2(3) GeV 2 . A similar analysis was performed for SU(3) in [11] and the con- 
clusions favour a D(0) = 0, although a finite nonvanishing D(0) was not completely excluded. The 
simulations use the Wilson action and different lattice spacings and volumes. For the SU(2) simu- 
lations the authors used a ~ 0.22 fm and volumes up to (27 fm) 4 , whereas the SU(3) simulations 
were performed with a ~ 0.10 fm and volumes up to (6.5 fm) 4 . 

Using the lattice data described before for j3 = 5.7, we are now in position of review the scaling 
analysis for the SU (3) gauge theory — see also [|l2|]. Note that for the largest volumes 64 4 - 96 4 , 
which are data from the Berlin-Moscow-Adelaide group, we only have access to D(0). 

Anyway, following JT0| , 11], we assume that in (4.1) the different functions depend on the 
lattice volume as A/V a . Then, it follows that an a > 1 means D(0) = in the infinite volume. The 
fits of the lattice data to the small set of volumes (L < 8.1 fm) give 



<M(0)> 2 


D(0)/V 


(M(0) 2 ) 


a 1.0537(50) 


1.0504(45) 


1.0530(50) 



and confirm the conclusions presented in jjllj], as we found a > 1. If one wants to use the full set 
of p = 5.7 lattices, we can only investigate the scaling behaviour of D(0)/V. In this case, the fit 
gives a = 1.0538(28) with a X 1 / d.o.f. = 0.87, which is in excelent agreement with the estimation 
using the smaller set of lattices, and again it suggests a D(0) = in the infinite volume. The 
j8 = 5.7 lattice data and the corresponding fits to A/V a can be seen in figure ^. Again, as in [O], 



if one assumes that the dependence with the volume of the different members of (4.1) are given 
by C/V + D/V a , then the data is well described by this functional form; in this sense, a finite and 
non- vanishing value for D(0) in the infinite volume is not excluded. 
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5. Power-law Behaviours From Ratios Of Propagators 

In [|13|] a method was proposed which tests the compatibility of the lattice data with a power 
law behaviour and, simultaneously, suppresses the finite volume effects. Since the method is not 
wellknown, for completeness, we will review it. The infrared propagator can be investigated using 
on-axis momenta which are defined as 

^N = -sin(^— J, « = 0,1,...,- (5.1) 

for a symmetry L 4 lattice. If the gluon propagator is described by a pure power law, i.e. 

D(q 2 ) =Z{q 2 f K -\ (5.2) 

one can define the following ratio 

R[n] 



The compatibility of the lattice 

R[n] = 2KR q [n] + C, (5.5) 



In 



q 2 [n+\}D{q 2 [n+1]) 
q 2 [n]D(q 2 \n}) 

~q 2 [n+\ 



2KR r ,\n] = 2k In 



q 2 [n] 



(5.3) 
(5.4) 



data with a power law can be tested fitting the lattice data to 
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Figure 8: Ratios of propagators for j3 = 5.7 lattices. 



where C is a constant which resumes both the deviations from a power law behaviour and the finite 



volume effects which are not eliminated by taking the ratio of propagators in R[n}. In [|13|], the 
authors analyzed such type of fits to the lattice data for large asymmetric 8 3 — 18 x 256 lattices 
— see figure [J. For the gluon propagator the fits give, within one standard deviation, K > 0.5 with 
K ~ 0.53 and a constant C which seems to approach zero as we go from 8 3 x 256 to 18 x 256. The 
first result suggests a D(0) = 0, while the second result suggests that C mainly resumes the finite 
volume effects. In what concerns the ghost propagator, figure [7] shows that either the data is still far 
from the linear behaviour or the ghost propagator does not follow a pure power law in the infrared 
region. 

The ratios for the /3 = 5 .7 propagators can be seen in figure |8[ We have tried to fit the data 
coming from the largest lattice volumes to (5.5) but it turns out that the % 2 /d.o.f '. were always 



too large, i.e. well above 2. This is probably due to rotational invariance violations as the infrared 
Berlin-Moscow-Adelaide data mixes different types of momenta which have different types of cor- 
rections due to the lack of rotational invariance. Indeed a similar effect is seen in results obtained 
from other lattices when one mixes on-axis with off-axis momenta. Given that we cannot distin- 
guish between the two types of points in the Berlin-Moscow-Adelaide data, we show in figure ^ the 
ratios for some asymmetric lattices and the large volume SU(2) propagators. The figure includes 



the K measured from fitting the lattice data to (5^). The 128 4 results for the SU(2) gauge group 
is the only data which gives a K < 0.5, within one standard deviation. However, the infrared 128 4 
propagator, see figure [K| shows larger fluctuations than any other calculation. Furthermore, it is 
one of the very few simulations where one sees an enhancing of the propagator in infrared region. 
Certainly, this is due to the small statistics. 



9 



Lattice Infrared Landau Gauge Gluon Propagator 



O. Oliveira 



1.5 



0.5 



- 

- 
- 

- 




, , 1 


- 


^^^^ K 






0.533(19) 


■ 16 3 x 256 






0.524(13) 


• 18 3 x 25 6 






0.539(18) 


♦ 80 4 [SU(2) P = 2.2] 






0.495(13) 


128 4 [SU(2)] 






, , 1 , , , , 1 , , 


, 1 , , , , 1 , 



O.d 



Figure 9: Ratios of propagators for SU(3) asymmetric lattices at j3 = 6.0 and SU(2) simulations. 



6. How Different are the SU(2) and SU(3) Infrared Gluon Propagators? 



As a final topic, we would like to discuss how different are the SU(2) and SU(3) propagators. 
The discussion of section |3| suggests that there should be a difference in the infrared. Indeed, the 
mass scales for SU(2) and SU(3) don't seem to be equal with M Sf/ ( 2 ) > ^su(3)- Moreover, the 
scaling analysis of the Cucchieri-Mendes bounds seems to give different conclusions when we use 
different gauge groups — see section |] and [10, 11]. However, in [14, 15] it was shown that the 



two propagators are similar for momenta larger than ~ 800 MeV. For smaller momenta, one can 
see some differences whose origin was not clear. If one uses now the SU(2) data from the S. Carlos 
group and the SU(3) Berlin-Moscow-Adelaide data one can compare results for volumes up to 
(17 fm) 4 . In order to compare the SU(2) and SU(3) propagators, they were renormalized by the 
condition 

D ^)lf-^ = J2 (6-D 

choosing /I = 3 GeV. The renormalization constants were computed after fitting D(q 2 ) to the 1- 
loop perturbative result for q > 2.5 GeV. In these fits for SU(3) we used the conic cut data, while 
for the SU(2) we used only the diagonal momenta. The renormalization constants were computed 
from the condition ( |dTT| ) using the results of the fits. The statistical error in the renormalization 
constants being around 10%. 

The renormalized propagators, which can be seen in figure [R], show clear differences in the 



infrared region. This result is not in contradiction with [J14J], which had a limited access to the low 
momenta region. Furthermore, reading D(0) from figure [H] it follows that 



D(0)su(2) j>_ = 2 
O(0)st/(3) ~ V.5 3 ' 



(6.2) 
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Renormalized Gluon Propagator -[1 = 3 GeV 




Figure 10: Comparing SU(2) and SU(3) gluon propagators. 



suggesting a D(0) su ^ ~ N c . This result is not necessarily in conflict with the large N c expansion 
because the gluon propagator is not a gauge invariant object. 

We would like to call the reader's attention that, although the propagators are different, the 
infrared exponents computed by the ratio method discussed in the previous section seem to be 
similar. 



7. Conclusions 

The results discussed in this article show that results obtained from lattice simulations for the 
gluon propagator in Landau gauge can be made compatible with both types of Schwinger-Dyson 
solutions, i.e. it can be seen either as a scaling like solution, where 0(0) = 0, or as a decoupling 
solution, where 0(0) is finite and non- vanishing. The raw lattice data shows a plateau at small 
momenta and, in this sense, it does look more like as a decoupling type of propagator. 

The analysis of the bounds derived in JTQ| ] are not conclusive when applied to the SU(3) sim- 
ulations [|ll|, [12|]. According to the analysis shown here for the /3 =5.7, the leading behaviour is 
already captured when one includes volumes as "small" as ~ (8 fm) 4 . Moreover, the performed 
scaling analysis is in good agreement with the ratios analysis, favouring a 0(0) = 0. This does not 
necessarily mean that the lattice simulations point towards a scaling like solution. Remember that 
the ratio method applied to the asymmetric lattice data shows that the ghost propagator does not 
follow a power law. In what concerns the gluon and ghost propagators, it looks like that before 
having a clear understanding of the finite lattice volume/spacing effects, it will be quite difficult to 
give a definitive answer on the nature of the propagators in the deep infrared region. 

Finally, the Landau gauge gluon SU(2) and SU(3) propagators are compared for momenta 
below 800 MeV. It turns out that there are clear differences and a scaling law D(0) su ^ ~ N c 
seems to hold. 
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